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Abstract 

The vacuum expectation values of the energy-momentum tensor are investigated for 
massless scalar fields satisfying Dicichlet or Neumann boundary conditions, and for the elec- 
tromagnetic field with perfect conductor boundary conditions on two infinite parallel plates 
moving by uniform proper acceleration through the Fulling-Rindler vacuum. The scalar 
case is considered for general values of the curvature coupling parameter and in an arbitrary 
number of spacetime dimension. The mode-summation method is used with combination of 
a variant of the generalized Abel-Plana formula. This allows to extract manifestly the con- 
tributions to the expectation values due to a single boundary. The vacuum forces acting on 
the boundaries are presented as a sum of the self-action and interaction terms. The first one 
contains well known surface divergences and needs a further regularization. The interaction 
forces between the plates are always attractive for both scalar and electromagnetic cases. An 
application to the 'Rindler wall' is discussed. 



PACS number(s): 03.70. +k, ll.lO.Kk 



1 Introduction 

The imposition of boundary conditions on a quantum field leads to the modification of the spec- 
trum for the zero-point fluctuations and results in the shift in the vacuum expectation values for 
physical quantities such as the energy density and stresses. In particular, vacuum forces arise act- 
ing on constraining boundaries. This is the familiar Casimir effect. The particular features of the 
resulting vacuum forces depend on the nature of the quantum field, the type of spacetime manifold 
and its dimensionality, the boundary geometries and the specific boundary conditions imposed on 
the field. Since the original work by Casimir in 1948 [1] many theoretical and experimental works 
have been done on this problem, including various types of boundary geometry and non-zero tem- 
perature effects (see, e.g., [2, 3, 4, 5, 6, 7, 8] and references therein). Many different approaches 
have been used: mode summation method with combination of the zeta function regularization 
technique. Green function formalism, multiple scattering expansions, heat-kernel series, etc. An 
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interesting topic in the investigations of the Casimir effect is the dependence of the vacuum char- 
acteristics on the type of the vacuum. It is well known that the uniqueness of vacuum state is 
lost when we work within the framework of quantum field theory in a general curved spacetime 
or in non-inertial frames. In particular, the use of general coordinate transformation in quantum 
field theory in fiat spacetime leads to an infinite number of unitary inequivalent representations 
of the commutation relations. Different inequivalent representations will in general give rise to 
different pictures with different physical implications, in particular to different vacuum states. 
For instance, the vacuum state for an uniformly accelerated observer, the PuUing-Rindler vacuum 
[9, 10, 11, 12], turns out to be inequivalent to that for an inertial observer, the familiar Minkowski 
vacuum. Quantum field theory in accelerated systems contains many of special features produced 
by a gravitational field avoiding some of the difficulties entailed by renormalization in a curved 
spacetime. In particular, near the canonical horizon in the gravitational field, a static spacetime 
may be regarded as a Rindler-hke spacetime. Note that, as it has been shown in Ref. [13], there is 
a class of solutions to the Einstein equations with a plane-symmetric matter distribution for which 
the corresponding external geometry is described by the Rindler metric ('Rindler walls'). Another 
motivation for the investigation of quantum effects in the Rindler space is related to the fact that 
this space is conformally related to the de Sitter space and to the Robertson- Walker space with 
negative spatial curvature. As a result the expectation values of the energy-momentum tensor 
for a conformally invariant field and for corresponding conformally transformed boundaries on 
the de Sitter and Robertson-Walker backgrounds can be derived from the corresponding Rindler 
counterpart by the standard transformation (see, for instance, [14]). 

In this paper we will consider the vacuum expectation values of the energy-momentum tensors 
for a scalar and electromagnetic fields in the region between two parallel plates moving by constant 
proper acceleration through the Fulling-Rindlcr vacuum. This problem for a single plate case 
was considered by Candelas and Deutsch [15] and by one of us [16]. In Ref. [15] the cases 
of conformally coupled Dirichlet and Neumann massless scalar and electromagnetic fields are 
investigated in the region of the right Rindler wedge on the right from the barrier. In Ref. [16] 
both regions, including the one between the barrier and Rindler horizon are considered for a 
massive scalar field with general curvature coupling parameter and Robin boundary conditions in 
arbitrary number of spacetime dimensions, and for the electromagnetic field. As in Ref. [16] (see 
also [17, 18, 19, 20, 21]), our regularization scheme here is based on a variant of the generalized 
Abel-Plana formula derived in Appendix A. This allows to extract form the vacuum expectation 
values the single boundary parts and to present the "interference" parts in terms of strongly 
convergent integrals useful for numerical evaluations. We have organized the paper as follows. In 
the next section we evaluate the vacuum expectation values of the energy-momentum tensor for the 
Dirichlet scalar. The corresponding interaction forces between the plates are investigated in section 
3. Section 4 is dedicated to the case of the Neumann boundary conditions. Then the vacuum 
densities and interaction forces for the electromagnetic field are considered in section 5. Section 
6 concludes the main rcsTilts of the paper and an application to the 'Rindler wall' is discussed. 
In Appendix B wc consider the case of the scalar field in two spacetime dimensions separately. 
An alternate representation of the vacuum expectation values for the energy-momentum tensor 
is obtained in Appendix C. 
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2 Vacuum energy-momentum tensor for a Dirichlet scalar 



We consider a real massless scalar (fi{x) field with general curvature coupling parameter ^ satisfying 
the field equation 

V^V^V' + (Rip = 0, (2.1) 

with R being the scalar curvature for a + 1-dimensional background spacetime, is the 
covariant derivative operator associated with the metric g^i,. For minimally and conformally 
coupled scalars C = C = {d ~ l)/4'i, respectively. By using field equation (2.1) it can be 
seen that the corresponding energy-momentum tensor (EMT) can be presented in the form 



(c - \) ^M.V,V^ - CV^V. - CR^, 



^\ (2.2) 



where Rfj,^ is the Ricci tensor. 

Let {ipa{x),^a{x)} is a complete set of positive and negative frequency solutions to the field 
equation (2.1), where a denotes a set of quantum numbers. Expanding field operator over these 
eigenfunctions and using the commutation relations it can be easily seen that the vacuum expec- 
tation values (VEV's) of the EMT are presented in the form 

{0\T,,\0) = ^T^,{ip^,cp*J, (2.3) 

a 

where for a scalar field the quadratic form T^i, {/, g} directly follows from the classical EMT given 
by Eq. (2.2). 

Our main interest in this paper will be the vacuum expectation values (VEV's) of the EMT 
in the Rindler spacetime induced by two parallel plates moving with uniform proper accelera- 
tion when the quantum field is prepared in the Fulling-Rindler vacuum. For this problem the 
background spacetime is fiat and in Eqs. (2.1), (2. 2) we have R — 0, R^i, — 0. As a result the 
eigenmodes are independent on the curvature coupling parameter and the EMT VEV's will de- 
pend on this parameter through the expression (2.2) only. In the following it will be convenient 
to introduce Rindler coordinates (r, ^,x) related to the Minkowski ones, (t,a;^,x) by 

t — ^ sinh T, — ^ cosh r, (2.4) 

where x = (x^, . . . ,x'^) denotes the set of coordinates parallel to the plates. In these coordinates 
the Minkowski line element takes the form 



ds^ = ^^dr^ - d^^ - d^\ (2.5) 

and a wordline defined by ^, x = const describes an observer with constant proper acceleration 

Assuming that the plates are situated in the right Rindler wedge x^ > \t\ we shall let the 
surfaces ^ = and C, = C,2, > Ci represent the trajectories of these boundaries, which therefore 
have proper accelerations ^ and (see Fig. 1). First we will consider the case of a scalar field 
satisfying Dirichlet boundary condition on the surface of the plates: 

1^=6= k=6= (2-6) 

To evaluate the VEV's of the EMT by Eq. (2.3) we need the form of the eigenfunctions ^Pai^). For 
the geometry under consideration the metric and boundary conditions are static and translational 
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Figure 1: The {x^,t) plane with the Rindler coordinates. The heavy hnes ^ = and ^ = ^2 
represent the trajectories of the plates. 



invariant in the hyperplane parallel to the plates. It follows from here that the corresponding part 
of the eigenfunctions has the standard plane wave structure: 



V^a = C(f){^) exp [i (kx - ujt)] , 



a 



{k2 



(2.7) 



The equation for 0(^) is obtained from field equation (2.1) on background of metric (2.5) and has 
the form 

eV"(o + e0'(o + - k\') m = 0, (2.8) 

where the prime denotes a differentiation with respect to the argument, and k = |k|. In the 
region between the plates the corresponding linearly independent solutions to equation (2.8) are 
the Bessel modified functions Uuik^) and Ki^{k^). The solution satisfying boundary condition 
(2.6) on the plate ^ = ^2 is in form 



DUk^, A;6) = I^Uk^2)KUkO " KUk^2)I^u.m. 



(2.9) 



Note that this function is real, Di^^lkC,, k^2) = D_i^{k^, k^2)- From the boundary condition on the 
plate ^ = ^1 we find that the possible values for uj are roots to the equation 



DUkii.ki2) = Q. 



:2.10) 



This equation has an infinite set of solutions. We will denote them by = uj£,ni i^on > 0, 
n = 1,2,..., and will assume that they are arranged in the ascending order uon < ^^Dn+i- 
The coefficient C in formula (2.7) is determined from the standard Klein-Gordon orthonormality 
condition for the eigenfunctions which for metric (2.5) takes the form 



dx / —(fa dr fa' 



5n 



(2.11) 



It can be easily seen that for any two solutions to equation (2.8), ^^^HO; 'm = 1,2 the following 
integration formula takes place 



(2.12) 



Taking into account boundary condition (2.6) from Eq. (2.11) for the normalization coefficient 
one finds 



C 



1 



Now substituting the eigenfunctions 

(p^{x) = CoDi^jy^ik^, k^2) exp [i (kx - ujonr)] 



(2.13) 



(2.14) 



into Eq. (2.3) and integrating over the directions of k for the VEV's of the EMT we obtain 
diagonal form (no summation over i) 



n=l huj{kC2) 

where \0d) is the amplitude for the Dirichlet vacuum between the plates, and 

1 



Ad 



(2.15) 



(2.16) 



In formula (2.15) for a given function G{z) we use the notations 



dG{z) 



dz 



dG{z) 



dz 

2 



+ ^-^\Giz)f + 

z dz 



1 



1 ^ 

- 2C + ^ 



+ 2C 



--4:IGWr + ^(i-^l|G(^)r, 



z dz 
2C-i 

^ 2 



\G{z)\', (2.17) 



(2.18) 



dG{z) 



dz 



+ I--. \G{z 



; i = 2,...,d,(2.19) 



where G{z) = Dii^{z, k^2), and the indices 0,1 correspond to the coordinates r, ^ respectively. It 
can be easily seen that for a conformally coupled scalar the EMT (2.15) is traceless. 

For the further evolution of VEV's (2.15) we will apply to the sum over n summation formula 
(A. 5) derived in Appendix A by making use of the generalized Abel-Plana formula [17]. This 
yields 



re 

dkk'^ / 

Jo 



duj 



sinh nuj 



TT 



Ukii)F^^[D^{ki,ki2)] 
Uk^2) D^{ki^M2) 



where we have introduced the notation 



(2.20) 



(2.21) 



and the functions F^'^^[G{z)], i — 0, 1, . . . , d are obtained from the functions f^'^^[G{z)] (see Eqs. 
(2.17)-(2.19)) replacing uj iu: 



F«[G(z)]=/«[G'(z), 



a; — > zo; . 



(2.22) 



The vacuum energy density, e, effective pressures in perpendicular, p, and parallel, p±, to the 
plates directions are determined by relations (no summation over i) 



e={On\T^\OD), p = -{Qd\T^\Od), P± = -{OdWIOd) , i = 2,...,d. 



(2.23) 
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It can be easily checked from Eqs. (2.20), (4.6) and (2.17)-(2.19) that they satisfy the standard 
continuity equation for the EMT, which for the geometry under consideration takes the form 

(2.24) 

For a conformally coupled scalar wc have an additional zero-trace relation e — p — {d — l)p± = 0. 
Let us consider the limit ,^2 ^ oo of general formula (2.20) for fixed ^. It can be easily seen that 
in this limit the VEV's take the form 

hm (0^|7;'=|0^) = (O^IT^'^IO^) + {Tt)^n''\^i,0, C > 6, (2-25) 

where 

(QjilThOR) ^ dkk'^ dcJsmhTTUjf'^Ki^ik^)] (2.26) 

TT Jo Jo 

are the corresponding VEV's for the FuUing-Rindler vacuum without boundaries, and the term 

{Tt)^D\Ci,0 = -AJ^ rdkk'^ r d^¥^,F^^iK^i^O] (2.27) 

Jo Jo K^[k^i) 

is induced in the region ^ > C,i hj the presence of a single plane boundary located at ,^ = ,^1. 
Expressions (2.27) are finite for all values ^ > ^-^id divergences are contained in the purely 
FuUing-Rindler part (2.26). These divergences can be regularized subtracting the corresponding 
VEV's for the Minkowskian vacuum. The subtracted VEV's 

{Tt)lul = {ORlTtlOn) - (Om^IOm) (2.28) 

are investigated in a large number of papers (see, for instance, [15, 16, 22, 23, 24, 25, 26, 27, 
28, 29, 30, 31, 32] and references therein). The most general case of a massive scalar field in an 
arbitrary number of spacetime dimensions has been considered in Ref. [28] for conformally and 
minimally coupled cases and in Ref. [16] for general values of the curvature coupling parameter 
(for the corresponding Green function see [22]). The formulae relevant to this paper are given in 
[16]. For a massless scalar VEV's (2.28) can be presented in the form 

V-'i /sub — orf-1_rf/2T-/J/o\ L „2-7raj , / T\d l^Z.Zy; 



2d-%<i/2r(rf/2) Jo e^'^'^ + (-1) 

(the expressions for the functions g^'^\uj) are given in Ref. [16]) correspond to the absence from 
the vacuum of thermal distribution with standard temperature T = (27r^)~^. As we see from Eq. 
(2.29), in general, the corresponding spectrum has non-Planckian form: the density of states factor 
is not proportional to uj'^~^duj. The spectrum takes the Planckian form for conformally coupled 
scalars in (i = 1, 2, 3 with g^^\ijS) = —dg'^'^\u) = 1, i = 1,2, . . . d,. It is interesting to note that for 
even values of spatial dimension the distribution is Fermi-Dirac type (sec also [33, 34]). For the 
massive scalar the energy spectrum is not strictly thermal and the corresponding quantities do 
not coincide with ones for the Minkowski thermal bath. 

The boundary induced quantities (2.27) are investigated in Ref. [15] for a conformally coupled 
d — 3 massless Dirichlet scalar and in Ref. [16] for a massive scalar with general curvature coupling 
and Robin boundary condition in an arbitrary number of dimensions. The single boundary part 
(2.27) diverges at the plate surface = ^^i with leading terms proportional to — for 
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i = 0,2, . . . ,d and to for i = 1 (see below). These leading terms vanish for a conformally 

coupled scalar, and for i = 0,2, ... ,d coincide with the corresponding quantities for a plane 
boundary in the Minkowski vacuum [16]. 

Now we turn to the limit ^ in formula (2.20), when the left plate coincides with the right 
Rindlcr horizon. In this limit in the second term on the right of formula (2.20) the subintegrand 
behaves as ^f"^ and tends to zero. As a result one obtains 



lim{OD\ThOD) = dkk'^ du smhnu f^\D,^{ki,ki2)]- (2.30) 

6^0 TT Jo Jo 



These quantities coincide with the corresponding ones induced in the region ^ < ^2 by a single 
plate at ,^ = ,^2- They are investigated in Ref. [16], where it has been shown that the VEV's (2.30) 
can be presented in the form similar to Eq. (2.25): 

lim(0^|7;'=|0^) = {Or\T^\Or) + {T^)^S'\i2,0, C < 6, (2.31) 
51^0 

where the expressions for the boundary part ^i)^D\'i2,C) the region ^ < ^2 are obtained from 
formulae (2.27) by replacing (see Ref. [16]) 

I^^K^, K^^I^, ^1^6, (2.32) 
By using Eqs. (2. 20), (2. 30), (2. 31) the parts in the VEV's induced by the existence of boundaries, 

{T^)^S = {MT^M - (Oi^mlO^), (2.33) 



can be written as 



Jo Jo Iuj[k^2) D^[k^i,k^2) 

In Appendix C we show that the VEV's (2.20) can be also presented in the form (C.4). 
Substituting Eq. (C.ll) into this formula, the boundary VEV's can be also written in the form 

Jo Jo K^{k^i) D^{k^i,kt,2) 

This expression is obtained from Eq. (2.34) by replacements (2.32). The case d = 1 needs a 
separate consideration and is investigated in Appendix B. It can be seen that the corresponding 
formulae for the VEV's are also obtained from the formulae given above in this section replacing 

POO 1 

Ad dkk'^-^^-, k^O. (2.36) 

Jo TT 

Now let us present the VEV's (2.20) in the form 

(oiT^'^lo),, = (Oii|T;1o^) + (7;'=)g''^(ei,0 + (^')S^'^(6,0+A(T;^)i.(6,6,0, 6 < e < 6, (2.37) 

where 

F«[D4A;e,A;6)] F«[if^(A;0]" 



A{T^')d^-AJ'1 dkk" dculUkCi 
Jo Jo 



UkC2)DUkCi,kC2) K^ikCi) 



(2.38) 
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is the 'interference' term. The surface divergences are contained in the single boundary parts and 
this term is finite for all values < ^ < ^2- An equivalent formula for A{T!^)d is obtained from 
Eq. (2.38) by replacements (2.32). In the hmit — > ^2 expressions (2.38) are divergent and for 
small values of ,^2/'Ci ~ 1 the main contribution comes from the large values of uj. Introducing a new 
integration variable x = k/00 and replacing Bessel modified functions by their uniform asymptotic 
expansions for large values of the order (see Ref. [35]) at the leading order over 1/(^2 — ^1) 
receives (no summation over i) 



for the single boundary terms, and 

1 



d{Cc - or (^) 



i = 0,2,...,d, 



1,2 



D 



d 



X 



f 

Jo 



dtf^ 



(c-Ce)(e2-ei 

22d-i^d/2r(rf/2) 



-d-l 



X 



(2.39) 
(2.40) 

(2.41) 



exp 



)r(^ 



+ exp t 



6-6 



D 



D, 



2,3,..., 



(2.42) 



for the 'interference' terms. Here Cr{s) is the Riemann zeta-function. Expressions (2.39), (2.41), 
(2.42) coincide with the corresponding formulae for two parallel plates geometry in d + 1 - di- 
mensional Minkowski spacetime with separation ^2 — 6 (see Ref. [36] for the conformally coupled 
case and Ref. [18] for the general case of the curvature coupling parameter Q. Note that in the 
limit under consideration the 'interference' term (2.42) for the vacuum perpendicular pressure 
dominates the single boundary induced terms, given by Eq. (2.40). 



3 Interaction forces between the plates 

Now we turn to the interaction forces between the plates. The vacuum force acting per unit 
surface of the plate at ^ = 6 is determined by the J-component of the vacuum EMT at this point. 
The corresponding effective pressures can be presented as a sum of two terms: 

pg^=pg\+pgU' ^ = 1'2- (3-1) 
The first term on the right is the pressure for a single plate at = 6 when the second plate is 
absent. This term is divergent due to the well known surface divergences in the subtracted VEV's. 
The second term on the right of Eq. (3.1), 

is the pressure induced by the presence of the second plate, and can be termed as an interaction 
force. For the plate at ^ = 6 the interaction term is due to the second summand on the right of 
Eq. (2.20). Substituting into this term ^ — ^2 and using the Wronskian relation for the modified 
Bessel functions one has 

(2) U C\ r AUUd-2 r A I^jk^l) /oQN 
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By a similar way from Eq. (2.35) for the interaction term on the plate aX ^ = we obtain 



P2(int)te'6) 



2^f Jo Jo 



duj 



(3.4) 



As the function D^{k^, k^2) is positive for < ^2, interaction forces per unit surface (3.3) and 
(3.4) are always attractive. They are finite for all ,^1 < ,^2, and do not depend on the curvature 
coupling parameter (. In the limit ,^1 C,2 these forces diverge due the contribution from the 
large values lv and in this limit by introducing a new integration variable we can replace the Bessel 
modified functions by their uniform asymptotic expansions for large values of the order. At the 
leading order for the perpendicular vacuum pressures we obtain formula (2.42) which corresponds 
to the standard Casimir attraction force for two parallel plates in Minkowski vacuum. 
From expressions (3.3) and (3.4) it follows that 



Pnlint) iCi , 6) > Pniint) (6 , 6) ■ 



(1) 



(3.5) 



This can be proved by using that the function z'^Ii^{z)K^{z) is monotonic increasing. The latter 
directly follows from the relations 



'1 + 



^2 



/ (a; + 1)2 1 

'1 + ^ — + - > 




(3.6) 
(3.7) 



The proof for the right inequalities in Eqs. (3. 6), (3. 7) is presented in Ref. [15]. The left inequalities 
are obtained from the recurrence relations for the Bessel modified functions. For instance, in the 
case of the function I^{z) one has: 




UJ 

+ - = 
z 



- + 



1 



00 

+ - > 

z 



+ 



'1 + 



(u; + ll 



(3.8) 



where we have used the right inequality in Eq. (3.6). The left inequality in Eq. (3.7) can be 
proved in a similar way. 

To see the monotonicity properties of functions (3.3) and (3.4) note that 



K\ TT^ — -?2- 



96 



96 



266 



/ dkk''-^ / 
Jo Jo 



duj 



Dl{kiiM2) 

.(1) f„(2) 
(int) \PD{mt)/ 



(3.9) 



It follows from here that for a fixed value of 6 (6) the quantity p^^lint) (PD(int)) monotonic 
increasing (decreasing) function on 6 (6)- By taking into account that both this quantities are 
negative we conclude that the modulus of the interaction force on the plate at ,^1 (,^2) is monotonic 
decreasing (increasing) function on 6 (6) ^ fixed value of 6 (6)- Ftom formula (3.3) it follows 
that 



96 



-{d + i)p 



D(int) 



96 



(3.10) 
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For i = 2 the both terms on the right are positive and hence, the same is the case for the function 
on the left. Therefore for a fixed the function p^jj^t) is monotonic increasing on ^2 and the 
modulus of the corresponding interaction force is monotonic decreasing function on ^2- In the 
case i = 1 the terms on the right in this formula have different signs. For a fixed value of ,^2 the 
function p^lj^t) is monotonic increasing on ,^1 near the horizon, ^1 — 0, and monotonic decreasing 
near the second plane, ^1 — > ^2- It follows from here the modulus of the corresponding interaction 
force acting on the plate at ^1 has minimum for some intermediate value. 

In the limit ^2 ^ Ci? introducing in Eq. (3.3) a new integration variable x — k^2, and making 
use the formula 



1 + 0{y' 



y = ^6/6, 



and the standard relation between the functions K^^ and one finds 



(2) 
Pd (mt) 



48e|+Mn2(26/ei) 



°° dxx'^ ^ 



x) 



1 + 



\nx 



(in(26/ei). 



The similar calculation for Eq. (3.4) yields 



Pd (int) 



24e2^-^e?ln^(2e2/ei) 



dxx Ko{x) 



1 + 



\ylx 



(in(26/ei 



(3.11) 



(3.12) 



(3.13) 



We have carried out numerical evaluations for the interaction forces by making use of formulae 
(3.3) and (3.4). In Fig. 2 the corresponding results are presented for C2^^PD(int)' ^ = 1,2 in the 
case d = 3 as functions on ^1/^2- 




0.1 0.2 0.3 0.4 0.5 0.6 



Figure 2: The d = 3 vacuum effective pressures determining the interaction forces between Dirich- 
let parallel plates, multiplied by ^tPolmt) (curve a) and C2PD(int) (curve b) as functions of the 
ratio {1/^2- 



4 VEV's and the interaction forces for the Neumann scalar 



In this section we will consider VEV's for the EMT in the case of a scalar field satisfying the 
Neumann boundary condition on the plates ^ = ^1,^2: 



dip I 



= 0. 



(4.1) 
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The corresponding scheme is similar to that given above for the Dirichlet case. The eigenfunctions 
to the field equation (2.1) have form (2.7) with 

As in the Dirichlet case this function is real. From the boundary condition on the plate ^ = we 
obtain that the corresponding eigenfrequencies are solutions to the equation 

NUkir, ^6) = lL{ki2)KUkii) - KUki,)Il{kii) = 0. (4.3) 

We will denote them by a; = ujjsim n = 1, 2, ... , arranged in the ascending order up^ri < ^Nn+i- The 
normalization coefficient C can be found from orthonormality condition (2.11) using integration 
formula (2.12): 

Substituting the eigenfunctions into the mode sum formula (2.3) one obtains 



roo °° I' (kf ) 

{On\T,'\On) = vrAA' / dkk'' E 77--T§^feM;) /^^^^-(^e, ^6)]U=..„, (4.5) 



where |Ojv) is the amplitude for the Neumann vacuum state between the plates, and the functions 
/*-*''[G(2;)] are defined in accordance with Eqs. (2.17)-(2.19). To sum the series over the eigen- 
frequencies u}]sin we will apply the summation formula derived in Appendix A, Eq. (A. 11). This 
yields 



- ^ r sinhTTo; ,,,,,,, Wi) F»[iV^(fce, k^. 



{Or,\Tt\0^) ^ A,5N dkk' dco r-^f'^[NUkC,kC2)] 

Jo JO TT 



2) 



[ TT ILik^2) Nl^{k^i,k^2) 

with functions F^'^^[G{z)] defined as in Eq. (2.22), and we use the notation 



(4.6) 



kC2) = KUkO - ^§4^W). (4.7) 

To identify the terms in Eq. (4.6) let us consider limiting cases. In the limit ^2 ^ oo, from Eq. 
(4.6) one obtains 

hm {0^\T,'\On) = (O^mlO^) + {Tt)^^'\^i,0, (4.8) 
52^00 

where the term 

{Tt)'^N\Ci, = -AaS^ r dkk' r du^^^F^\K^{ki)] (4.9) 

Jo Jo Kl[k^i) 

is induced in the region ,^ > ,^1 by a single Neumann boundary located at ,^ = ^1. This quantities 
for d = ?) case are investigated in Ref. [15]. In the limit ^1 the left plate coincides with 
the Rindler horizon and the second term in the figure braces in Eq. (4.6) vanishes. In this case 
the VEV's coincide with the corresponding expressions for a single plate at ^ = ^2 induced in the 
region ^ < ^2- They are investigated in Ref. [16], where it has been shown that the VEV's (2.30) 
can be presented in the form similar to Eq. (4.8): 

lim(0^|7;>^) = {Qr\T^\Qr) + {T^)f\i2,i), e < 6, (4.10) 
41^0 
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where the expressions for the boundary part {Tj<^i'^)N^\^2,0 region ^ < ^2 sue obtained 

from formulae (4.9) by replacements (2.32). 

By using Eqs. (4. 6), (4. 10) the parts in the VEV's induced by the existence of boundaries, 

{Tt)^^^ = {ON\Tt\0^) - {OR\Tt\OR), (4.11) 

can be presented as 

ffi')i5' (6. &. = Ut)T (&. - r ak r ^M^A^^f^m^. (4.12) 

Similar to the Dirichlet case, the Neumann boundary VEV's can be also written in the form 

(T^fs (6.6.0 ^ (i;'>r(6.6 - A^st fikk' rd.§§4^»|M, (4.13) 

Jo Jo K'^[kt,i) N^{k^i,K2) 

with {T^)n^\^i,^) being the VEV's induced by a single Neumann boundary located at ^ = ^1. 
As we see, this expression is obtained from (4.12) by replacements (2.32). 
Now let us present the VEV's (4.6) in the form 

(o^|7;^|o^) = (o^|7;^|o«)+(7;'=)S'')(6,0+(^')S'^(6,0+A(7;')iv(6,6,0, 6 < C < 6, (4.i4) 

where 



A{T^')r,^-AJ^ dkk^ dcuIiikCi) 
Jo Jo 



im2)K{k^^M2) K{k^,] 



(4.15) 



is the 'interference' term. An equivalent formula for A{TI')n is obtained from Eq. (4.15) by 
replacements (2.32). 

'Interference' term (4.15) is finite for all < ^ < ^2, ^1 < '^2? and diverges in the limit ^1 — >■ ^2- 
In this limit the main contribution into the cu-integral comes from the large values u. Introducing 
a new integration variable x — k^\/ui and using the uniform asymptotic expansions for the Bessel 
modified functions in the leading order one obtains that the quantities A{T^)n coincide with the 
VEV's for two parallel plates in c? + 1-dimensional Minkowski spacetime with separation ,^2 — 
[36, 18]. The corresponding expressions are given by formulae (2. 41), (2. 42) with the opposite sign 
of the integral term on the right of formula (2.41). 

Now we turn to the Neumann vacuum effective pressures determining the forces acting on 
the plate due to the presence of the second plate (interaction forces). This force acting per unit 
surface of the plate ^ = ^21 Pyvlint) defined by the J-component of the second term on the right 
of formula (4.12) at = .^2- The nonzero contribution comes from the last term on the right 
of Eq.(2.18) (with replacement (2.22)). Using the standard Wronskian relation for the Bessel 
modified functions one obtains 

.gu(6,6)^^f«"r-fii±^. 

By a similar way for the interaction force per unit surface of the first plate from the second term 
on the right of Eq. (4.13) at ^ = ^1 we receive 

PnU^^.^-) - 2eJo io K{ki.)K{ki,M2) ■ ^'-'^^ 
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Note that pressures (4. 16), (4. 17) are independent on the curvature couphng parameter. It can be 
seen that the function I'^{z) / K'^{z) is monotonic decreasing, and as a result N^{k^i,k^2) < for 
^1 < In combination with Eqs. (4.16), (4.17) it follows from here that pjv*(int) < 0, i = 1, 2, and 
hence, as in the Dirichlet case, the Neumann interaction forces are always attractive. By using 
that the function z'^I'^{z)K'^{z) / {z"^ + uj^) is monotonic decreasing (this can be proved by using 
inequalities (3. 6), (3. 7) ) we see that 

pSU(^i'^2)>p!iU(6,6)- (4.18) 

In the limit ^2 replacing the Bessel modified functions by their uniform asymptotic expansions 
we can see that to the leading order from Eqs. (4. 16), (4. 17) the standard Casimir interaction force 
is obtained for two parallel plates with separation ^2 — ^1 the d + 1-dimensional Minkowski 
spacctime. 

From formulae (4.16), (4.17) one has 

= = ^ f f ^(i±^im^^±^. ,4,,) 

9^2 0^1 2^i5 Jo Jo K {k^i, k^2) 

As seen from here for a fixed value of ,^2 {^1) the modulus of the interaction force acting on the 
plate at = ^2 (^1) is a monotonic increasing (decreasing) function on ^1 (,^2)- For the other partial 
derivatives, similar to the Dirichlet case, one has the relation (3.10) with replacement D ^ N. 
In particular, we can see that dp\^^^^^^^/ d^2 > 0. The Neumann effective pressures determining 
the interaction forces per unit surface given by Eqs. (4.16), (4.17) are plotted in Fig. 3 as 
functions of ^1/^2 for the case 0? = 3. As seen from Fig. 2 and Fig. 3 the Dirichlet and Neumann 




0.1 0.2 0.3 0.4 0.5 0.6 



Figure 3: The d = ?> vacuum effective pressures determining the interaction forces per unit surface 
between Neumann parallel plates, multiplied by ^2Pi)(int) (curve a) and ClPi^jint) (curve b) as 
functions of the ratio ^1/^2- 

vacuum interaction forces are numerically close to each other. This is a consequence of that the 
subintegrands in formulae (3.3) and (4.16) and in formulae (3.4) and (4.17) are numerically close. 
This can be also seen analytically by using relations (3. 6), (3. 7). 
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5 Electromagnetic field 



We now turn to the case of the electromagnetic field in the region < ^ < ^2- We will assume 
that the mirrors are perfect conductors with the standard boundary conditions of vanishing of 
the normal component of the magnetic field and the tangential components of the electric field, 
evaluated at the local inertial frame in which the conductors are instantaneously at rest. By 
considerations similar to those given in Ref. [15] for d — 3, it can be seen that the corresponding 
eigenfunctions for the vector potential may be resolved into one transverse magnetic (TM) 
and d — 2 transverse electric (TE) (with respect to ^-direction) modes A^^, a = 0,1, ... ,d — 2, 
a = (k, a;): 

;-e9/9e,-We,0,...0)^0a, ^ = 1, TMmode, (5.1) 

(5.2) 



AfJ' 



0, 2, . . . , - 2, TE modes. 



where the polarization vectors obey the following relations 



el = el = 0, e.^, = -fcXa', e^A;^ = 0. (5.3) 

From the perfect conductor boundary conditions one has the following conditions for the scalar 
fields (^o-qI 



</7(Ta|^=6 ^ VTo\i=ii =0, (7 = 0, 2, . . . , d - 2, 



«=6 = 0. 



(5.4) 



As a result the TE/TM modes correspond to the Dirichlet/Neumann scalars. In the correspond- 
ing expressions for the eigenfunctions A^^ the normalization coefficient is determined from the 
orthonormality relation 

rfx t = -— 5aa'(^aa'. (5.5) 

On the base of this normalization condition for the separate scalar modes one has 

27rV2 

(/'era = -^Cz'Zi^^^{ki, ki-i) exp \i (kx - lozut)] , (5.6) 

where Z — D iox a — 0,2, . . . ,d — 2 and Z — N ior a — 1, and the coefficients Co and Cat are 
defined in accordance with Eqs. (2. 13), (4.4). Substituting the eigenfunctions (5.1), (5.2) into the 
mode sum formula 



mm = E I'dkY^T^A^^^A^^^ 



a=0 



(5.7) 



with the standard bilinear form for the electromagnetic field EMT one finds 



r(¥) 



dkk' J2 P'rT.Cy^^[Ziu.zM,k^2)], Po^d-2, A = l, (5.8) 



o-=0,l 



n=l 



where (3o and f3i are the numbers of the independent polarization states for TE and TM modes 
respectively. In Eq. (5.8) for a given function G{z) the following notations are introduced 



&G{z)] = 



dG{z) 



dz 
dG{z 



1 + ^ \G{z) 



dz 



d-5 
d-1 



\G{z)\' + 



z^ 

d-3 
d-1 



\G{z) 
dG{z) 



(5.9) 



dz 



-^\G{z) 



2,3,. ..,d. 
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By making use of the summation formulae derived in the Appendix A the VEV's are presented 
in the form 

(5.10) 

where I^^ — /j^, /^^^ = and the same notations for the functions K^^, Z^^. The functions F^W 
are obtained from Eqs. (5.9) replacing cu — > iuj: 

FB[G{z)]^fi^[G{z),u;^iu;]. (5.11) 

It can be easily checked that the components (5.10) obey the covariant conservation equation and 
the corresponding EMT is tracclcss for d = 3. The first term in the figure braces of Eq. (5.10) 
corresponds to the VEV induced by a single plate at ^ = .^2 in the region < ^2- For the case 
d = 3 they arc investigated in Rcf. [16]. The generalization for an arbitrary d is straightforward 
and these quantities are presented in the form 

mW\i2.i) = {MT^\^r) - -AAa / dkk' duY. (3^ Z)),T F^[Uka (5.12) 

where (0r|7^'^|0r) are the VEV's for the Fulling-Rindler electromagnetic vacuum without bound- 
aries. By the way similar to that given in Ref. [16] for the case of a scalar field, it can be seen 
that 



(0«|r, [0^) - (OmIT, [Om) - 2d-%./2r(d/2) L e^- + (-l)-^ \\ 



'd- 1 -2/' ^ 



2tu 



+ 1 



, (5.13) 



where Im — d/2 — 1 for even d> 2 and l„i = (c? — l)/2 for odd d > 1, and the value for the product 
over I is equal to 1 for 0? = 1, 2, 3. In Eq. (5.13) we have introduced notations 

/ol(a;) = -rf/£l(a;) = 1 + (5.14) 
fit) (, ,\ _ f (1) (, ,\ , {d-'i-){d-3) 

JoemK^) — /0emi<^j H , l — Z,...,a. 

For physically most important case d = 3, formula (5.13) leads to the standard result derived by 
Candelas and Deutsch in Ref. [15]. 

An alternative form for the vacuum EMT in the region between two plates is 

(o|7;^|o) = (o^|7;^io^) + (o|7;^io)(^'')(6,0- 

- btA. r dkk^ r d. E (5.15) 

where (OlT^^lO)^^^^^, is the vacuum EMT induced by a single boundary at ,^ = ,^1 in the region 
^ > ,^1. The latter is obtained from (5.12) by replacements (2.32). For the interaction force p^^^^i^xy 
i — 1,2 per unit area of the plate at ^ = from Eqs. (5.10) and (5.15) one obtains 

- rdk¥-' Hdu, T (-ire ^fi."(fc6)(i + ^Vt'^;r ,5 

e,.., ^ -|ir--=r- E i-^rp.!^'^:,^. (.17) 

^^2-^0 Jo ^^01 Id,'{k^2) Z^'{k^i,k^2) 
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Recalling that (-1)'"^^'") > we see the electromagnetic interaction forces are attractive. Note 

that Pe'^(int) = {d- '^)p%^t) +PN{int)- I" hmit ^ and to the leading order over 1/(^2 - 
^i) from these expressions the electromagnetic Casimir interaction force between plates in the 
Minkowski spacetime is obtained. 

6 Conclusion 

It is well known that the uniqueness of vacuum state is lost when wc work within the framework 
of quantum field theory in a general curved spacetime or in non-inertial frames. In this paper 
we have considered vacuum expectation values of the energy-momentum tensor for scalar and 
electromagnetic fields between two infinite parallel plates moving by uniform proper acceleration, 
assuming that the fields are prepared in the Fulling-Rindler vacuum state. As the boundaries 
are static in the Rindler coordinates no Rindlcr quanta arc created and the only effect of the 
imposition of boundary conditions on quantum fields is the vacuum polarization. For the scalar 
case the both Dirichlet and Neumann boundary conditions are investigated. The VEV's are 
presented in the form of mode sums involving series over zeros cu — cuon or cu — cuNn oi the 
functions Dii^{k^i, k^2) and N-^{k$^i, k^2) respectively. To sum these series we derive in Appendix 
A summation formulae for these types of series using the generalized Abel-Plana formula. The 
application of these formulae allows to extract from the VEV's the parts due to the single plate. 
The latters are investigated previously in Refs. [15, 16]. The boundary induced parts are presented 
in two alternative forms, Eqs. (2.34), (2.35), for the Dirichlet case, and Eqs. (4. 11), (4. 12) for the 
Neumann case. Various limiting cases are studied. In particular, in the limit when the left plate 
coincides with the Rindlcr horizon the corresponding VEV's are the same as for a single plate 
geometry. The vacuum forces acting on boundaries contain two terms. The first ones are the 
forces acting on a single boundary then the second boundary is absent. Due to the well-known 
surface divergences in the VEV's of the energy-momentum tensor these forces are infinite and need 
an additional regularization. The another terms in the vacuum forces are finite and are induced 
by the presence of the second boundary and correspond to the interaction forces between the 
plates. These forces per unit surface do not depend on the curvature coupling parameter and 
are determined by formulae (3. 9), (3. 10) for the Dirichlet scalar and by formulae (4. 16), (4. 17) for 
the Neumann scalar, and are always attractive for both plates. In particular, they are the same for 
conformally and minimally coupled scalars. For given ^i, ^2 the modulus of the interaction force 
is larger for the plate at ^ = ^1 (see inequalities (3.5) and (4.18)). For small distances between the 
plates at the leading order the standard Casimir result on background of the Minkowski vacuum is 
rederived. The case of the electromagnetic field is considered with the perfect conductor boundary 
conditions in the local inertial frame in which the boundaries are instantaneously at rest. The 
corresponding eigenmodes are superposition of TE and TM modes with Dirichlet and Neumann 
boundary conditions respectively. The VEV's of the electromagnetic EMT in the region between 
the plates are given by formulae (5.10) and (5.15). The corresponding vacuum interaction forces 
per unit surface, Eqs. (5. 16), (5. 17), are obtained by summing the Dirichlet and Neumann d — ?> 
scalar forces, and are attractive for all values of the proper accelerations for the plates. 

The results obtained in this paper can be applied to the geometry of two parallel plates near 
the 'Rindler wall'. With the x coordinate perpendicular to the wall and with the {x^^x^) plane 
located at the centre of the wall, a; = 0, the static plane-symmetric line element can be written as 

ds"" = e'^^'^Ue - dx"" - e^(^)dx2, X = (x^ x^), (6.1) 
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where i>{x) and X{x) are even functions. For this metric the Einstein equations with the diagonal 
matter energy-momentum tensor T^'^'''^ — (iiag(£:(™), — — — admit two classes of 
solutions. For the first one A'(0) > 0, and the corresponding external solution (the solution in the 
region x > Xg, where T^^^^ = 0, with x = Xg being the boundary of the wall) is described by the 
standard Taub metric [37]. For the second class of internal solutions A'(0) < 0, and the external 
solution is presented by the metric 

dsl^^ = e"' [1 + 2nas{x - Xs)]'^dt'^ - dx^ - e^^dx^ (6.2) 

where Vs = ^^(2:^5), = A(a;s), and 

(7, = 2e-^^/2-^^ j^' + 3p("^)) e'^/^+^dx (6.3) 

is the mass per unit surface of the wall. For a given equation of state p*^™^ = p('")(£(™)) the 
parameters Vg^Xs^Xg are functions of the central pressure p'^™')|^=0) ^-nd ^"^^ determined by the 
internal solution of the Einstein equations (see Ref. [13] for the case of the equation of state 
corresponding to the incompressible liquid). Now redefining 

i{x) ^x-Xs + ^— , T = 27r(7,e''^/^i, e^'''^x' ^ x\ 2, 3, (6.4) 

2710-5 

from Eq. (6.2) we obtain the Rindlcr metric in the form (2.5). Hence, the VEV's for the EMT in 
the region between two plates located at a; = and x = X2, Xi > Xg near the 'Rindler wall' are 
obtained from the results given above substituting = i{xi), i = 1,2 and ^ = i{x). Note that for 
as > 0, X > Xg one has ^{x) > Ci^s) > and the Rindler metric is regular everywhere. 
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A Summation formulae over zeros of Diz and N^^ 

In this section we will derive a summation formula over zeros z — cuok oi the function 

Diz{x, y) = Ii^{y)Ki^{x) - Ii^{x)Ki^{y), y > x. (A.l) 

As we saw in section 2 the VEV's of the EMT for the Dirichlet scalar between two plates in the 
Fulling- Rindler vacuum are expressed in the form of series over these zeros. To derive a summation 
formula we use the generalized Abel-Plana formula [17]. Let us choose in this formula 

f{z) = —sinhnz F{z), (A.2) 

71 

{y)I-iz{x) + Iiz{x)I-iz{y) 

0^'^ = l^J^) 
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with a meromorphic function F[z) having poles z = Zk ii^ the right half-plane Re^; > 0. The sum 
and difference of functions (A. 2) are presented in the form 



g[z) ± f[z) = ^ , ^ . F{z). 



(A.3) 



By taking into account that the zeros ujr>k are simple poles of the function g{z) for the function 
R[f{z),g{z)\ in the generalized Abel-Plana formula one obtains 



R[f{z),g{z)] = 2m 



E 



F{z) 



\Z=UJDk 



+ 



(A.4) 



+ Ef^es 



Z=Zk 



Diz{x,y) 



-^isgn(Imzj,)z (y)-^-isgn(Imzfc)2 {x) 



where the zeros udu are arranged in ascending order. As a result we obtain the following summa- 
tion formula 



I-iz{y)Iiz{x) 
k=l 'd^-^iz{x,y) 



E 



Fiz) 



\Z=LODk 



1 /""^ 

— / svahixz Fiz)dz 

TT^ Jo 

F{z) 



(A.5) 



- E^-es 



Z=Zk 



F>iz{x,y) 



isgn{lrazk)z 



isgn{lm.zu)z 



[x) - 



— — / dz- — , - , , — Li x)I-z[y). 



,{y)K,{x) - I^,{x)KM 

Here the condition for the function F(z) is easily obtained from the corresponding condition in 
the Generalized Abcl-Plana formula by using the asymptotic formulae for the Bessel modified 
function and has the form 



\F{z)\ < e{\z\)e- 



y\ 2|Imz| 

xj 



Rez > 0, \z\ 



oo, 



(A.6) 



where |2;|e(|^|) —>■ when \z\ — >• oo. 

A similar formula can be obtained for the series over zeros z = cuNk, k = 1,2, . . . of the function 

Niix,y) = nMKix) - KMlUx), y > x. (A.7) 

For this let us substitute in the Generalized Abel-Plana formula [17] 

2i 



■sinh TTZ F(z) 



(A.8) 



^^^-j ^ i'iz{y)i'-iz{.x) + riz{.x)r_iz{y) ^^^^ 



NL{x,y) 

Using these expressions it can bcc easily seen that 



9{^) ± m 



2I',,{x)I'^M 



NL{x,y) 



Fiz). 



For the function R[f{z),g{z)] now one obtains 



R[m,9{z)] = 2ni 



,{y)I'iz{x) 



y- 



F{z) 



\Z=UlNk 



+ 



(A.9) 



(A.IO) 



+ ER-es 

k 



F{z) 



Z=Zk ATI 



N'iz{x,y) 



^isgn(Im.Zk)z{,y) —isga(hnzk)z{,x) 
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As a result wc obtain the following summation formula 

7:V , mi=... = i sinhnz F{z)dz- (A.ll) 



F{z) 

where the corresponding condition for the function F{z) has the form (A. 6). 



B d = 1 case: Direct evoluation 



For d = 1 case the linearly independent solutions to equation (2.8) are e^*'^^"^. The normalized 
eigenfunctions satisfying Dirichlet boundary conditions (2.6) are in form 



D 



where we use the notation 



sm an 



7TO 



a 



Tin 
ln(6/Ci) 



vrln(6/0 



(B.l) 



(B.2) 



ln(6/6) ■ 

Substituting eigenfunctions (B.l) into mode-sum formula (2.3) and applying to the sum over n 
the Abel-Plana summation formula one finds 

(ODmlO^) - (OMmlOM) = + (J;^)g''^(6,0 + C ^'i1''2^.;.^ diag(l,0) - 



27re^ 



TT 



c 



.i2in2(6/6) ' Hi2li\ 

Here the subtracted purely FuUing-Rindler part without boundaries, (T'i'^)(^b)' P^^^ ™" 

duced by a single boundary at ^ = ^2 are given by formulae [16] 



'27reW(6/e) 
+ 1 /-^ /-.X (oicot a — 1) 



diag(l,-l)(B.3) 



(^')(s1.l) = ^(C-l^)diag(l,-l), 



C 



-diag(l + l/ln(e/6),-l)- 



(B.4) 



(B.5) 



Note that the expression (B.5) for a single boundary part is valid for both regions { < and 
^ > ^2- Now for the vacuum interaction forces between the plates one obtains 



TT 



PDiint) 24efin2(e2/ei)' 



1,2. 



(B.6) 



In the limit d ^2 to the leading order we recover the standard Casimir result on background 
of the 2D Minkowski spacetime. 
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For the case of the Neumann boundary conditions (4.1) the normahzed eigenf unctions have 
the form 

(p^ ^ —= cos( an) , n = 0,1,2,..., (B.7) 

where u and a are given by the same relations (B.l) and (B.2) as in the Dirichlet case. The 
substitution of these cigcnfunctions into the mode-sum formula shows that the VEV's of the 
EMT for the Neumann boundary conditions can be obtained from the corresponding formula for 
the Dirichlet case, Eq. (B.3), replacing in the boundary part ( 



C Alternative representation for the VEV's 

As a solution to equation (2.8) satisfying first boundary condition (2.6) one could take the function 

DUk^, fcCi) = hUk^i)K,^{kO - K,^{ki^)lUkO- (CI) 

Now from the boundary condition on the plate = .^2 (2.6) wc find the possible values for u being 
roots to the equation (2.10). For the normalization coefficient we receive 

(27r)'^-^7..(fe6) ^^-£"^^-^ '^^"""' ^ ^ ^ 

The VEV's of the energy - momentum tensor are obtained in a diagonal form 

(0^|7;>^)=7r^,5f r^^^^f — (C.3) 



For the further evoluation of VEV's (C.3) we can apply to the sum over n summation formula 
(A.5). This gives 



Jo Jo I ^ 

I.^{k^2)F(^^[D^{kC,k^i)]\ 



(C.4) 



I-.ik^i)D^ik^i,k^2) J 
This form of the VEV's is equivalent to Eq. (2.20). To see this let us consider the quantities 

a« - - sinhTT^ fWfD- (k£ k£ )] - P i^- ^^^^^^^^-^^^^^ ^^^'^^ fC 5) 

5, -^J^ dujsmhn^f [D..{k^,k^,)\ /.(^0)/-.(^e,) ' ^ ^ ^ 

where j = 1,2. Two representations (2.20) and (C.4) will be equivalent if 

= (C.6) 

To prove this let us consider the difference 

<l2 -<A = - smh7ru;sW _ / duj S"- >, C.7 
n Jo Jo D^{k^i,k^2) 
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where we have introduced the notations 



s« = /«[A.,(fce,A:6)]-/^^^[A.(fce,^6)], (C.8) 



By using the standard relation between the Bessel modified functions it can be seen that the first 
integral in formula (C.7) can be presented as 

2 Jo h.{ki2)KUkii) - h.{kii)KUki2) ' ^'^ 

where the function s*^*^ iIiioik^2) Ki^{k^i) — Iiio{k^i) Ki^{k^2)) has no poles. For the term with the 
first (second) summand in the numerator rotating the integration contour by angle — 7r/2 (7r/2) 
in LV complex plane and noting that the integrals over arcs with large radius vanish (subintegrand 
behaves as we see that 



^ hUk^i)I-Uk^2) - I-.Uk^i)hUkC. 



2) Jo D^{k^,,k^2) ^ ' 



2 Jo lUki2)KUkii) - lUki^)KUki2 

Hence, the difference (C.7) is equal to zero, which directly proves Eq. (C.6) 

By taking into account Eq. (C.5) from Eq. (C.6) in the limit ^2 — 00 one obtains the following 
useful relation 

-/ da; sinh7ru;/W[Aa;(A;C,^Ci)] = -/ duj sinhiru f^'\Ki^{k^)]+ (C.ll) 

^Jo '''"\l-.{kir)K^{ki,) K^k^y ^""-^"^^i 
Substituting Eq. (C.ll) into Eq. (C.4) one finds formula (2.35). 
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